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I. INTRODUCTION 

An important part of the visual system of vertebrate animals are the neural connections between the eye and the brain. 
At an initial stage of ontogenesis the ganglion cells of the retina have random synaptic contacts with the tectum, a part of 
the brain which plays an important role in processing optical information. In the adult animal, however, neighboring retinal 
cells project onto neighboring cells of the tectum (see Figure [TJ. Further examples of these so-called retinotopic projections 
are established between the retina and the corpus geniculatum laterale as well as the visual cortex, respectively [1]. This 
conservation of neighborhood relations is also realized in many other neural connections between different cell sheets. For 
instance, the formation of ordered projections between the mechanical receptors in the skin and the somatosensorial cortex 
is called somatotopy. An even more abstract topological projection arises when the spatially resolved detection of similar 
frequencies in the ear are projected onto neighboring cells of the auditorial cortex. A further notable neural map in the auditory 
system was discovered in the brain of the owl, where neighboring cells of the Nucleus mesencephalicus lateralis dorsalis 
(MLD) are excited by neighboring space areas, i.e. every space point is represented by a small zone of the MLD [2]. The variety 
of examples suggest that there must be some underlying general mechanism for rearranging the initially disordered synaptic 
contacts into topological projections. 

In the early 1940s, Sperry performed a series of pioneering experiments in the visual system of frogs and goldfish UHH]. Fish 
and amphibians can regenerate axonal tracts in their central nervous system, in contrast to mammals, birds and reptiles. Sperry 
crushed the optical nerve and found that retinal axons reestablished the previous retinotopically ordered pattern of connections 
in the tectum. Then in the early 1960s Sperry presented his chemoaffinity hypothesis which proposed that the retinotectal map 
is set up on the basis of chemical markers carried by the cells @]. However, experiments over several decades have shown that 
the formation of retinotectal maps cannot be explained by this gradient matching alone [6]. 

The group of von der Malsburg suggested that these ontogenetic processes result from self-organization. The basic notion in 
their theory is the following: Once a fibre has already grown from the retina to the tectum, the fibre moves along by strengthening 
its contacts in some parts of its ramification and by weakening them in others. It is assumed that these modifications are 
governed by two contradictory rules (71 Is] : on the one hand, synaptic contacts on neighboring tectal cells stemming from fibres 
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FIG. 1: In the course of ontogenesis originally disordered mappings between retina and tectum evolve into ordered projections. 

of the same retinal region support each other to be strengthened. On the other hand, the contacts starting from one retinal cell 
or ending at one tectal cell compete with each other. In the case that retina and tectum are treated as one-dimensional discrete 
cell arrays, extensive computer simulations have shown that a system based on these ideas of cooperativity and competition 
establishes, indeed, retinotopy as the final configuration fl. This finding was confirmed by a detailed analytical treatment of 
Haussler and von der Malsburg [9] where the self-organized formation of the synaptic connections between retina and tectum 
is described by an appropriate system of ordinary differential equations. Applying the methods of synergetics II 1 OL ll 111 for 
one-dimensional discrete cell arrays, they succeeded in classifying the possible retinotopic projections and to discuss the criteria 
which determine their emergence. The more complicated case of continuously distributed cells on a spherical shell was partially 
discussed in Ref. Il2tl . 

It is the purpose of this paper to follow the outline of Ref. Il3tl and generalize the original approach by elaborating a model 
for the self-organized formation of retinotopic projections which is independent of the special geometry and dimension of the 
cell sheets. There are three essential reasons which motivate this more general approach. First, neurons usually do not establish 
1 -dimensional arrays but 2- or 3-dimensional networks. Hence the 1 -dimensional model of Haussler and von der Malsburg can 
only serve as a simplistic approximation of the real situation. Secondly, we want to include cell sheets of different extent, which 
is a more realistic assumption than neural sheets with the same number of cells. The third reason is that a general model is able 
to reveal what is generic, i.e. what is independent of the special geometry of the problem. Thus, here we generalize the Haussler 
equations to continuous manifolds of arbitrary geometry. By doing so, we proceed in a phenomenological manner and relegate 
a microscopic derivation of the underlying equations to future research. 

It should be emphasized that our main objective is not the biological modelling of retinotopy. Instead of that our consider- 
ations are devoted to the analysis of the dynamics of the nonlinear Haussler equations by using mathematical methods from 
nonlinear dynamics and synergetics. For the more biological aspects of retinotopy and the vast progress in modelling various 
retinotopic ally ordered projections during the last twenty years we refer the reader to the reviews Wl llil [l5ll . 

In Section [TT] we present the general framework of our model and introduce the equations of evolution for the connection 
weights between retina and tectum. We then perform in Section [III] a linear stability analysis for the equations of evolution 
around the stationary uniform state and discuss under which circumstances an instability arises. In Section [IV] we apply the 
methods of synergetics, and elaborate within a nonlinear analysis that the adiabatic elimination of the fast evolving degrees 
of freedom leads to effective equations of evolution for the slow evolving order parameters. They approximately describe the 
dynamics near the instability where an increase of the uniform growth rate of new synapses onto the tectum beyond a critical 
value converts an initially disordered mapping into a retinotopic projection. Finally, Section [V] and [yjprovide a summary and 
an outlook. 



II. GENERAL MODEL 

In this section we summarize the basic assumptions of our general model. 
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A. Manifolds and Their Properties 

We start with representing retina (R) and tectum (T) by general manifolds M.t and Ai r, respectively. In the framework of 
an embedding of these manifolds in an Euclidean space of dimension D, the coordinates Xr, xt of the corresponding cells can 
be represented by 

xr = (x R ,x 2 R , . . . ,x R ) , xr€Mr; xt = (x^,x T , • • • > xt & Mr ■ (1) 

In the following we need measures of distance, i.e. metrics g^ v , g^ v on the manifolds. The intrinsic coordinates of the d- 
dimensional manifolds Mr, M.t are denoted by r M , Thus, the vectors (fl~|i of the Euclidean embedding space can be 
parametrized according to xr = Xfl(r M ), xt = Xt^). With the covariant metric tensors 

r _ 9xr dx R rp _ dx T dx T 

the line elements on the manifolds are given by (cIsr) 2 — g^ v dr^dr v , (c?st) 2 = gj iv dt^dt 1 ' . The geodetic distances between 
two points of the manifolds read 



s R , = 



g% dr^dr- , sf t , = / JgT <W*df . (3) 



V 



We define a measure for the magnitudes of the manifolds by 

M T = Jdt, Mr = Jdr, 



(4) 



where we integrate over all elements of M.t, -Mr. We characterize the neural connectivity within each manifold M.t, -Mr 
by cooperativity functions cr(t, 1/), CR(r,r'). In lack of any theory for the cooperativity functions we regard them as time- 
independent, given properties of the manifolds which are only limited by certain global plausible constraints. We assume that 
the cooperativity functions are positive 

c T (t,t')>0, c R (r,r')>0, (5) 
that they are symmetric with respect to their arguments 

CT(t,t') = cr(t',t), c R {r,r') = c R (r \r) , (6) 
and that they fulfill the normalization conditions 

dt'c T {t,t') = 1, [dr'c R (r,r') = l. (7) 



Furthermore, it is neurophysiologically reasonable to assume that the cooperativity functions Or(t, t'), cr(t, r') are larger when 
the distance between the points t, t' and r, r' is smaller. This condition of monotonically decreasing cooperativity functions can 
be written as 

c T (t,t') > c T (t,t") if (4<) 2 < (4») 2 , c R {r,r') > c R (r,r") if (s*,) 2 < (s*,,) 2 . (8) 

B. Equations of Evolution 

The neural connections between retina and tectum are described by a connection weight w(t, r) for every ordered pair (t, r) 
with t 6 Mt, t € M.R. In this paper we are interested in the temporal evolution of the connection weight w(t, r) which 
is essentially determined by the given cooperativity functions or(t,t'), cr{t, r') of the manifolds A4t, -Mr. To this end we 
generalize a former ansatz of Haussler and von der Malsburg [9] and assume that the evolution is governed by the following 
system of ordinary differential equations 01611 : 



w(t,r) = a + w{t,r) jdt' \ dr' 'c T {t,t') c R {r,r') w{t' >') 
w(t, r) 



2M T 
w(t, r) 
' 2M R 



dt' 
dr' 



a + w(t',r) J dt" J dr'or(t',t")c R (r,r')w(t",r') 
a + w{t,r') jdt' fdr"cT(t,t')cR(r',r")w(t',r") 



(9) 



4 




FIG. 2: Illustrations for the respective contributions to the generalized Haussler equations i ll It . Discussion see text. 

Here a denotes the uniform growth-rate of new synapses onto the tectum, which will be the control parameter of our system. 
These equations of evolution represent a balance between different cooperating and competing processes. To see this, we define 
the growth rate between the cells at r and t 

f(t, r, w)=a + w(t, r) J dt' J dr'c T (t, t') c R (r, r') w{t' , r') , (10) 

so that the generalized Haussler equations (O reduce to 

^t,r) = f(t,r,w)-^£- Jdt'f(t',r,w)-^£ J ' dr> f(t,r',w) . (11) 

The cooperative contribution of the connection between r' and t' to the growth rate between r and t is given by the product 
w(t, r)cr(t, t')cji(r, r')w(t' , r') as shown in Figure^. Therefore, this cooperative contribution is integrated with respect to r', 
t' and added to the uniform growth rate a to yield the total growth rate (TTOb between r and t. Apart from this cooperative term in 
the equations of evolution (fTTt . the remaining terms describe competitive processes. The second term accounts for the fact that 
growth rates between r and t' compete with the connections between r and t (see Figure^). Correspondingly, the third term 
describes the competition of the growth rates between r' and t with the connections between r and t (see Figure^). 



C. Lower Limits for the Connection Strength 



Now we show that the evolution of the system due to the generalized Haussler equations (O leads to a lower bound for the 
connection weight. To this end we assume the inequality 

0<w{t,r)<W (12) 

to be fulfilled for some initial configuration. Then we conclude that the quantity 

C(t, r,w)= { dt' ( dr'c T (t, t') c R {r, r') w(t',r') (13) 
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is positive as both the cooperativity functions CT(t, t'), cr(t", r') and the connection weight wit' , r') are positive due to (0I and 
(fl2T >. On the other hand we read off from the normalization of the cooperativity functions (Q that C(t, r, w) cannot be larger 
than W: < C(t, r, w) < W. With this we can find a lower bound for w(t, t) as follows. The growth rate (TTOb reads together 
with©: f(t,r,w) = a + w(t, r) C(t, r, w). It can be minimized by setting C(t, r, w) = 0, i.e. 

f(t, r, U>)min = a, (14) 

whereas its maximum value follows from C{t, r, w) — W: 

fit,r,w) m ^ = a + W 2 . (15) 

To obtain a lower bound for w(t,r) in the Haussler equations ( fTTT ), we insert the minimum ( TT4| ) of the growth rate for the 
cooperative first term and its maximum dl~5T > for the remaining competitive terms: 

w(i,r) min = a- w{t,r) (a + W 2 ) . (16) 

Hence a small but positive wit, r) is prevented by a positive rate a from becoming zero. In this way we can conclude that the 
connection weight wit, r) is positive, when the inequality ( fT2b is valid in an initial configuration. All further investigations will 
concentrate on solutions of the Haussler equations (O with wit, r) > 0. Note that, in particular, the growth rates ( TTOb for such 
configurations are positive. 



D. Complete Orthonormal System 

To perform both a linear and a nonlinear analysis of the underlying Haussler equations © we need a complete orthonormal 
system for both manifolds M.t and M.r. With the help of the contravariant components <? T M , g R ^ of the metric introduced in 
Section lHAl we define the respective Laplace-Beltrami operators on the manifolds 

A T = -L d x (g^V^d^j , A R = -^d x (<$\/^) , (17) 



where gr, g R represent the determinants of the covariant components gj^, of the metric. The Laplace-Beltrami operators 
allow to introduce a complete orthonormal system by their eigenfunctions ip\ T it), ip\ R (r) according to 

At i>x T it) = xl T <Aa t it) , A R ^ Xr (r) = X x R ^ M • (18) 

Here At, denote discrete or continuous numbers which parameterize the eigenvalues x^ T , X\ R °f tne Laplace-Beltrami 
operators which could be degenerate. By construction, they fulfill the orthonormality relations 



dt^ XT it)ry T (t) = S XT y T , J dr^ XR ir)r X ' R (r) = 8 XRX>R , (19) 
and the completeness relations 

E>Xx (*') = *(*-*')> E^H( r )^ R (0 = <5(r-r'). (20) 

Note that the explicit form ( fTTb of the Laplace-Beltrami operators enforces the eigenvalues Xa t =o = > X\ R =o = with the 
constant eigenfunctions 

^ T = o(t) = 7W ' ^ R=o{r) = vfe (21) 

because of (|4|i and the orthonormality relations ( fT9l . The cooperativity functions can be expanded in terms of the eigenfunctions 
according to 

ct(M') -EE F XT y T 4, XT itWy T it') , c R (rS) = E E F ^'J^Wx R ir') . (22) 

At X' t Xr X' r 

In the following we assume for the sake of simplicity that the corresponding expansion coefficients are diagonal F\ T \i = 
f\ T S\ T \' T , F\ R \' R = f\J\ R \' R , so we have 

ct(M') = E h T ^ T mi T (t') , c R (r,r') = E /A H V>A fl (r)^ fl (r') ■ ( 2 3) 
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Thus, ip\ T (t), i/j\ R (r) are not only eigenfunctions of the Laplace-Beltrami operators as in dT8b but also eigenfunctions of the 
cooperativity functions according to 

dt' c T {t,t')ip\ T (t') = f\ T ip\ T {t) , J dr c R (r,r')ipx R (r') = f\ R ip\ R {r) . (24) 

Note that the normalization of the cooperativity functions (0 and the orthonormalization relations ( fT9b lead to the constraints 

/a t =o = f\ R =a = 1 • 

III. LINEAR STABILITY ANALYSIS 

Now we employ the methods of synergetics ITol [Till and investigate the underlying equations of evolution (O in the vicinity 
of the stationary uniform solution. Inserting the ansatz w(t, r) — wo into the Haussler equations ©, we take into account © 
as well as the normalization of the cooperativity functions (0. By doing so, we deduce wq = 1. Let us introduce the deviation 
from this stationary uniform solution v(t, r) = w(t, r) — 1 , and rewrite the Haussler equations ||9). Defining the linear operators 

C(t,r,x) = Jdt 1 Jdr' c T {t,t r ) c R (r,r') x{t' ,r') , (25) 

B(t,r,x) = — l — [dt'x(t',r) + —— fdr'x(t,r'), (26) 



2M T J 2M R 
the resulting equations of evolution assume the form 

v(t, r) = L(t, r, v) + Q(t, r, v) + K(t, r, v) . (27) 

Here the linear, quadratic, and cubic terms, respectively, are given by 

L(t,r,v) = -av + C(t,r,v)-B(t,r,v)-B(t,r,C(t,r,v)), (28) 
Q(t, r, v) = v (C(t, r, v) - B(t, r, v) - B(t, r, C(t, r, v))\ - B (t, r, v C(t, r, v)j , (29) 

K(t,r,v) = -vB (t,r,vC(t,r,v)^ . (30) 



To analyze the stability of the stationary uniform solution we neglect for the time being the nonlinear terms in (I27l i and investigate 
the linear problem 

v(t,r) = L{t,r,v). (31) 

Solutions of (|3~TT l depend exponentially on the time r, v(t, r) — v\ T \ R (t, r) exp (A\ T x R t) with v\ T \ R and A\ T \ R denoting the 
eigenfunctions and eigenvalues of the linear operator L: 

L(t,r,vx T \ R ) = Ax T \ R V\ T \ R (t,r) . (32) 

Now we use the complete and orthonormal system on the manifolds A4t, -Mr, which have been defined in Section iHDl and 
show that the eigenfunctions of L are products of the form 

v\ T x R (t,r)=ipx T (t)ipx R (r). (33) 

Indeed, when the operator (EBT l acts on (l33l . the expansion of the cooperativity functions ( f23l leads, together with the orthonor- 
mality relations (fT9l l. to 

C(t, r, vx T x R ) = h T As vx T x R (t, r) . (34) 
Thus, the operator C has the eigenfunctions vx T x R {t, r) with the eigenvalues fx T fx R ■ In a similar way we obtain for the operator 

v x t Xr At = A fl = , 

Uil.r. r VrXl: ) ={ V Xt xJ2 \ T = 0,\ R ^0;\ R = 0,\ T ^0, (35) 

otherwise . 



7 



Im A 

A 



stable modes Re A s <C 



/ 



unstable modes Re A u = 



-►Re A 



FIG. 3: Schematic representation of the eigenvalues J37b at the instability. The unstable part consists of those eigenvalues which nearly vanish 
whereas the stable part lies in a region separated by a finite distance from the stable part. 



Combining the eigenvalue problems d34l ), ( f33T > for C and B, we find 

ffx T fx R vx T x R \t = \r = , 

fx T fx R v XT x R /2 X T =0,X R ^0;X R = 0,X T ^0, (36) 
otherwise . 

Thus, we conclude from (I34il-(l3"6il that the linear operator L fulfills the eigenvalue problem (l32l l with the eigenfunctions (l33l 
and the eigenvalues 

{—a — 1 At = A;? = , 

-a + {fx T fx R -l)/2 \ T =0,\ RJ L0;\ R = 0,\ TJ L0, (37) 
-a + fx T fx R otherwise . 

By changing the uniform growth rate a in a suitable way, the real parts of some eigenvalues (|37| > become positive and the 
system can be driven to the neighborhood of an instability. Which eigenvalues (|37] l become unstable in general depends on 
the respective values of the given expansion coefficients fx T , fx R - The situation simplifies, however, if we follow Ref. H 
and assume that the absolute values of the expansion coefficients f\ T , fx R are equal or smaller than the normalization value 
fa = 1: \fx T | < 1 j \ fx R | < 1. Then the eigenvalue in d37l ) with the largest real part is given by some parameters A r , X R with 
A m ax = A A « A « = — a + fx u fx u - Thus, the linear stability analysis reveals that the instability arises at the critical uniform 
growth rate 

a c = Re(fx^fx R ) (38) 

and that its neighborhood is characterized by 

Re(A A¥A .)«0; Re(A ArA J«0, (A T ; X R ) ? (A T ; A^) . (39) 

Consequently, the absolute values of the eigenvalues of the unstable modes (A^; A^) are much smaller than those of the stable 
modes (A T ; X r ) ^ (A T ;A%): 

|Re (A A « A - )| « |Re (A AtA J| , (A T ; Xr) ? (A T ; A^) . (40) 
The resulting spectrum is schematically illustrated in Figure [3] 



IV. NONLINEAR ANALYSIS 

In this section we perform a detailed nonlinear analysis of the Haussler equations (O. Using the methods of synergetics ITol 
[Till we derive our main result in form of the order parameter equations which describe the emergence of retinotopic projections 
from initially undifferentiated mappings. 
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A. Unstable and Stable Modes 



We return to the nonlinear equations of evolution (|27] > for the deviation from the stationary uniform solution v(t, r). As the 
eigenfunctions ip\ T (t), i/j\ R (r) of the Laplace-Beltrami operators At, Ar represent a complete orthonormal system on the 
manifolds M.t, -Mr, we can expand the deviation from the stationary solution according to 

v(t,r)= V XTXR ^ XT (t)^ XR (r) . (41) 

Here we have introduced Einstein's sum convention, i.e. repeated indices are implicitly summed over. The sum convention is 
adopted throughout. Motivated by the linear stability analysis of the preceding section, we decompose the expansion (PiTt near 
the instability which is characterized by I 



v(t,r) = U(t,r) + S(t,r). (42) 

We can expand the unstable modes in the form 

U(t,r) = U X u Xl ^(t)^ X u(r), (43) 

where the expansion amplitudes f7 A « a% will later represent the order parameters indicating the emergence of an instability. 
Correspondingly, 

S(t,r) = S\ T \ R ipx T (t)^>X R (r) (44) 

denotes the contribution of the stable modes. Note that the summation in (l44l is performed over all parameters (At; \r) except 
for (A T ; A^), i.e. from now on the parameters (At; Xr) stand for the stable modes alone. In the following we aim at deriving 
separate equations of evolution for the amplitudes J/a^A" > S XtXr . To this end we define the operators 



(x) := JdtJ dr ^ (r) x(t, r) , (45) 

Kx R (x) := [dt [ dr rx T (tWx R W r) , (A T ; \r) ^ (A T ; X U R ) , (46) 



which project, out of v(t, r), the amplitudes of the unstable and stable modes, respectively: U\uyi = P A « a" (^) ; S XtXr — 
P\ T \ R (v) . These equations follow from (l42b — (f46b by taking into account the orthonormality relations (fT9] l. With these projectors 
the nonlinear equations of evolution (ETT i decompose into 

Ux^ T x R = A A?A . U X u A « + P X u Xl (Q(t, r,U + 5)) + P X u X u (k{t, r,U + S)) , (47) 

S XtXr = Ax t x r S Xt x r + Px t Xr (Q(t, r,U + 5)) + P Xt x r (k{t, r,U + 5)) . (48) 

Note that we used the eigenvalue problem (l32i l for the linear operator L and its eigenfunctions (l33l to derive the first term on 
the right-hand side in d47| i and ( |48l >, where Einstein's sum convention is not applied. 

In general, it appears impossible to determine a solution for the coupled amplitude equations d47l i. (US). Near the instability 
which is characterized by d38l . however, the methods of synergetics ITolfTlll allow elaborating an approximate solution which 
is based on the inequality d40b - To this end we interpret (|4Qb in terms of a time-scale hierarchy, i.e. the stable modes evolve on a 
faster time-scale than the unstable modes: 

1 1 

» T s - ^ ■ (49) 



|Re(A A . A .)r |Re(A ATA J| 

Due to this time-scale hierarchy the stable modes S XtXr quasi-instantaneously take values which are prescribed by the unstable 
modes U x ^x R - This is the content of the well-known slaving principle of synergetics: the stable modes are enslaved by the 
unstable modes. In our context it states mathematically that the dynamics of the stable modes S XtXr is determined by the center 
manifold H according to 

Sx T x R = H XtXr \Ux^x r ) ■ (50) 

Inserting (l50l in (|48| | leads to an implicit equation for the center manifold H which we approximately solve in the vicinity of 
the instability below. By doing so, we adiabatically eliminate the stable modes from the relevant dynamics. Then we use the 
center manifold H in the equations of evolution d47| i, i.e. we reduce the original high-dimensional system to a low-dimensional 
one for the order parameters Ux u x u ■ The resulting order parameter equations describe the dynamics near the instability where 
an increase of the uniform growth rate a beyond its critical value d38l converts disordered mappings into retinotopic projections. 
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B. Integrals 

It turns out that the derivation of the order parameter equations contains integrals over products of eigenfunctions which have 
the form 



7 a(i)a(2)...a(") = / dxtpl{x)i/j xl t)(x)tp x ^)(x) ••• ^ A (n)(x), (51) 
where A, x stand for the respective quantities At, t and Xr, r of the manifolds M.t and Mr. Examples for such integrals are: 
I x = Jdx ift(x) , = Jdx V A (z) i>x> (as) , 4 a» = /<& V> A (z) <Aa' (x) V>a» (x) . (52) 
The first two integrals of (l52l follow from the orthonormality relations (TL9l l by taking into account d2"H : 

/ A = VM(5ao , 4 = 5 xx' , (53) 

where M corresponds to Mt or Mr, respectively. Note that we will later make frequently use of the following consequence of 
d52> and d53j: 

fdxi)i(x) = 0. (54) 

Integrals with products of more than two eigenfunctions cannot be evaluated in general, they have to be determined for each 
manifold separately. At present we can only make the following conclusion. Expanding the product ipy (x) tpy (x) in terms of 
the complete orthonormal system 

ipx> (x) tp\" (x) = Cx>x"x"'ip\>" (x) , (55) 

the latter integral of ( 1521 is given by 



'A' A 

In addition, we will need also integrals of the type 



I x'x" = Ca'A"a • (56) 



• / A(DA(2)...A(") = / dx tp XW (x) ll) X (_ 2 ) (x) ■ ■ ■ V>A(») (x) , (57) 

for instance, 

Jaa' = J dx4>\(x) ip\>(x) . (58) 
Again we use the orthonormality relations (fT9l l. the expansion (|53T >, and take into account (fJTJ to obtain 

J AA < = \^Cavo, (59) 
where again M corresponds to Mt or Mr, respectively. 

C. Center Manifold 

Now we approximately determine the center manifold (T50b in lowest order. To this end we read off from ( |29l , d30i l. and d48b 
that the nonlinear terms in the equations of evolution for the stable modes S\ T \ a are of quadratic order in the unstable modes 
Ux^x^- Thus, the stable modes can be approximately determined from 

Sx T x n = A AtAr Sx T x n + Nx T x n (U) (60) 

with the nonlinearity 

Nx T x R (U) = Px T x R (UC(U) - UB(U) - UB (c(U)) - B (uC(U)) ) . (61) 
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Using the definitions of the linear operators (l25l l. d26b and the decomposition of the unstable modes (T43b as well as the projector 
for the stable modes d46l ), we see that the second and the third term in doTT i vanish due to (l54l 



P\ T \ R (UB(U)) =Px T x R (UB(C(U))) =0 



whereas the first term yields 



and the fourth term leads to 



P: 



\ tXr (uc(u)) = fx™,' h«. ill xy u x *x R u K , x ~, 



Px T x R {B(UC{U)\ \ = -/a-A-E/a-a-CVa- 



Therefore, we read off from (f6Tb— (f64b the decomposition 



N x T x R {U) - Qa|a|,A«'A^' Ux«\»Uxu>x%' , 



where the expansion coefficients are given by 



QXtXr r r 

X"X R ,X^'X R > - /AS/JA*' 



^AS AS,' ^aI A" ' - 77 ( nrr- ^A£A£' ^A«A^' ^A T 



1 / 1 



\ t a t R R 2 VvMr 



+ vfc Ja « a * Sxr0 ) 



(62) 



(63) 



(64) 



(65) 



(66) 



Note that Einstein's sum convention is not to be applied. To solve the approximate equations of evolution for the stable modes 
(l60l with the quadratic nonlinearity in the order parameters (l65l >. we assume that the center manifold (l50l l has the same quadratic 
nonlinearity: 



Sx T x R - H\Zx*,\«>\»> Ux^x R Uxyx R ' ■ (67) 
Inserting (j67j in (1601 . we only need the linear term in ( |47l i to determine the expansion coefficients of the center manifold: 



H 



A«A£,A£'A&' ~ (^ a t a h + ^X^'X R ' ^X T X R ) ^ x™X R ,X™' X R 



(68) 



Here, again, Einstein's sum convention is not to be applied. Therefore, the Eqs. 
of the center manifold. 



define the lowest order approximation 



D. Order Parameter Equations 

Knowing that the center manifold depends in lowest order quadratically on the unstable modes near the instability, we can 
determine the order parameter equations up to the cubic nonlinearity. Because of (l29l i. d30l l. and d47l i they read 



Ux^x R — Aa«a% Ux™x% + Nx*x a (U, S) , 
where the nonlinear term decomposes into three contributions: 

N x , Xr (U, S) = Qa»a« (U) + Kr^x^U) + K 2 ^ X ,(U, S) 



(69) 



(70) 



The first and the second term represent a quadratic and a cubic nonlinearity which is generated by the order parameters them- 
selves 



Qx%x H {U) = Px«x r (uc(u)-ub(u)-ub(c(u))-b(uc(u))) , 

Ki,x>x~(U) = ~Px^ t x r (ub(uC(U))) , 



(71) 
(72) 
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whereas the third one denotes a cubic nonlinearity which is affected by the enslaved stables modes according to 

K 2 , KX u(U,S) = P KX u(uc(s)-ub(s)-ub(c(s))-b(uc(s)) 



+SC{U) - SB(U) - SB C(U) )-B[ SC(U) 



(73) 



It remains to evaluate the respective contributions by using the definitions of the linear operators d251 l. (f26b and the decomposi- 
tions (1431 . d44l as well as the projector (05]). We start by noting that the last three terms in dTTT l vanish due to d54b . i.e. 



Px^xi (UB(U)) = P x * Xr (i/B (C(U))) = P X u Xl (B (UC(U))) =0, 
so the first term in ( FTTb leads to the nonvanishing result 

Qx^x R {U) = f x «» / A «» I x Z, X u„ ^a|'a«" u ^'^r' u K">- u r" ■ 
Correspondingly, we obtain for (l72l ) 



Ki,x^ R (U) 



2 f XV^'" f A^ 



U X «' Xr ' U X yx R " U\^>i>\v 



\» ' \" "A" "' ^A" A" ' J\%, "A" " 

A H A K A H T T T T 



Furthermore, taking into account d54l ), we observe that four of the eight terms in d73l vanish: 

Pa«a S (SB(U)) ,P X u X u (B (UC(S))) ,P X u Xl (SB (c(U))) ,Pa«a* (p (SC(C7))) =0. 
The nonvanishing terms in d73l read 

Px£A£ = fx T fx R lly\ T 7 Af'A n U X u, X u,S XtXr , 



P. 



'a«a% ySC(U)j - f X y I^, At ^a|'a r u K'>~r' s x t Xr 



and 



P: 



•a«a S (c^5(5) 



2 VVW 
1 



<5a t #A«A«' ^aC'Ar 



as well as 



/V S (i/s(c(S))) 



•i (-L 
1 



5x T S X u X u, f Xn i A |/ As 



+ 77jfT £\ r q£a*A r ' /At I \l'X T ) U X™'X R > Sx T X R 



It 



where we used /o = 1 in the last equation. Therefore, we obtain for ( l73l 



K. 



,X«X R (U, S) - U X ui Xr > S Xt X r { [fx T fx R + /A£' /a^' ] ^'At 7 Af'A R 



l r l 

2 



■S Xt0 S X u X u, (l + fx R ) I^ Xr + 



<Ja r o%a S ' (1 + /a t ) /a«'A 3 



'Mr 



(74) 



(75) 



(76) 
(77) 

(78) 
(79) 



(80) 



(81) 



(82) 



Taking into account d67| ), we read off from (1691 . d75l ), ( fTot , and d82l that the general form of the order parameter equations is 
independent of the geometry of the problem: 
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few mode amplitudes U\u 

of the slow 
linear unstable modes v\u 







order parameter 
slaving principle 
center manifold S = h(U) 



order parameter equations 

+N X "(U,h(U),a) 



i 


many mode amplitudes S A s 
of the fast 
linear stable modes v\* 







FIG. 4: Circular causality chain of synergetics for the order parameter equations of the generalized Haussler equations (5). The control 
parameter a denotes the growth rate of new synapses onto the tectum. 



The corresponding coefficients can be expressed in terms of the expansion coefficients f\ T , f\ R of the cooperativity functions 
d23l l and integrals over products of the eigenfunctions ip\ T (t), ip\ R (r) which have the form (fSTb or 1571 . They read 



and 



X R X R 'X R " 



JAJ" JX R " ^A£'A£" 1 \^ t '\ R > 



(84) 



Mr 



1 jX R r 



x Jx^'xi'" j + < [fx T fx R + fx%> fx%> ] Ixyx T I X\'Xr 



An 



Sx T oSx^' (1 + fXn) 'A S 



+ &*0 ^A R A«' (1 + At) Ixj'Ai 



At A_r 



(85) 



As is common in synergetics, the coefficients (1851 1 in general consist of two parts, one stemming from the order parameters 
themselves and the other representing the influence of the center manifold H. 



With (I83ll-(l85ll we have derived the generic form of the order parameter equations for the connection weights between two 
manifolds of different geometry and dimension. These equations represent the central new result of our synergetic analysis. 
Specifying the geometry means inserting the corresponding eigenfunctions of the Laplace-Beltrami operators dTvT > into the 
integrals (IBTT l. d57l > appearing in (l84l and (l85l l. Because the synergetic formalism needs not be applied to every geometry anew, 
our general procedure means a significant facilitation and tremendous progress as compared to the special approach in Ref. 



V. SUMMARY 

In this paper we have proposed that the self-organized formation of retinotopic projections between manifolds of different 
geometries and dimensions is governed by a system of ordinary differential equations (0 which generalizes a former ansatz 
by Haussler and von der Malsburg Jgt]. The linear stability analysis determines the instability where an increase of the 
uniform growth rate a beyond the critical value d38l converts an initially disordered mapping into a retinotopic projection. 
Furthermore, it gives rise to a decomposition of the deviation from the stationary uniform solution v(t, r) near the instability in 
unstable and stable contributions. By inserting this decomposition in the nonlinear Haussler equations ||9), we obtain equations 
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for the mode amplitudes of the unstable and stable modes, respectively. In the vicinity of the instability point the system 
generates a time-scale hierarchy, i.e. the stable modes evolve on a faster time-scale than the unstable modes. This leads to 
the slaving principle of synergetics: the stable modes are enslaved by the unstable modes. In the literature this enslaving 
5 = h(U) is usually achieved by invoking an adiabatic elimination of the stable modes, which amounts to solving the equation 
5 = 0. However, the mathematically correct approach for determining the center manifold h(U) is to determine it from the 
corresponding evolution equations for the stable modes lfl7ll . It can be shown that only for real eigenvalues this approach leads 
to the same result obtained by the approximation 5 = 0. Thus, it is possible to reduce the original high-dimensional system 
to a low-dimensional one which only contains the unstable amplitudes. The general form of the resulting order parameter 
equations (l83l is independent of the geometry of the problem. It contains typically a linear, a quadratic and a cubic term of 
the order parameters. As a general feature of synergetics, the coefficients (|83l , (l85l l consist of two parts, one stemming from 
the order parameters themselves and the other representing the influence of the center manifold on the order parameter dynamics. 

Our results can be interpreted as an example for the validity of the circular causality chain of synergetics, which is illustrated 
in Figure @] On the one hand, the order parameters, i.e. the few amplitudes of the slowly evolving linear unstable modes 
v\u , enslave the dynamics of the many stable mode amplitudes 5a» of the fast evolving stable modes v\s through the center 
manifold. On the other hand, the center manifold of the stable amplitudes acts back on the order parameter equations. 

VI. OUTLOOK 

The order parameter equations d83l-d85ll represent the central new result of this paper, and in the forthcoming publication 
lfl8ll they will serve as the starting point to analyze in detail the self-organization in cell arrays of different geometries. To this 
end we assume that the manifolds are characterized by spatial homogeneity and isotropy, i.e. neither a point nor a direction 
is preferred to another, respectively. This additional assumption requires the manifolds to have a constant curvature and their 
metric turns out to be the stationary Robertson-Walker metric of general relativity [ 19]. We therefore have to discuss the three 
different cases where the curvature of the manifolds is positive, vanishes, or is negative. This corresponds to modelling retina 
and tectum by the sphere, the plane, or the pseudosphere. 

A further intriguing problem concerns the question under what circumstances non-retinotopic modes become unstable and 
destroy the retinotopic order. One could imagine that some types of pathological development in animals corresponds to this case. 

As already mentioned, lacking any theory for the cooperativity functions, we have regarded them as time-independent given 
prop erties of the manifolds. They are determined by the lateral connections between the cells of retina and tectum, respectively 
02011 . But neither a reason for their time-independence nor a detailed discussion of their precise mathematical form is available. 
To fill this gap it will be necessary to elaborate a self-consistent theory of the cooperativity functions. 

Our generalized Haussler equations are fully deterministic. In real systems, however, there are always fluctuations. To take 
into account such unpredictable small variations a stochastic force has to be added to the deterministic part of the equation. 
Such fluctuations are known to play an important role, especially in the vicinity of instability points 1I2TI l22tl . 

Finally, delayed processes could be included in our considerations. Synergetic concepts have been successfully applied to 
time-delayed dynamical systems in Refs. |[T7i I23I l24l l25i l26ll . In neurophysiological systems delays occur due to the finite 

j 1 

propagation velocity of nerve signals H27L 12811 as well as the finite duration of physiological processes such as the change of 
synaptic connection weights. Thus, it would be also worthwhile to expand the investigations to time-delayed Haussler equations. 
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